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L°°-ERROR ESTIMATE FOR THE FINITE ELEMENT METHOD ON 
TWO DIMENSIONAL SURFACES 


HEIKO KRONER* 

Abstract. We approximate the solution of the equation 
(0.1) -Asu + u = f 

on a two-dimensional, embedded, orientable, closed surface S where —As denotes the Laplace Bel¬ 
trami operator on S by using continuous, piecewise linear finite elements on a triangulation of S with 
flat triangles. We show that the L°°-error is of order 0{h^ \ logh|) as in the corresponding situation 
in an Euclidean setting. 
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1. Introduction. During the last years several articles appeared which deal 
with the numerical solution of linear partial differential equations which are defined 
on a hypersurface in Roughly spoken their common aim is to show that concepts 
and properties which are well-known in an Euclidean setting carry over to the surface 
case. Without claiming completeness we summarize some steps towards this goal. 

In [B] the finite element approximation of the Laplace-Beltrami equation 

(1.1) -Asu = f 

on a surface S with continuous, piecewise linear elements (on a polyhedral approxi¬ 
mation of S) is presented and it is shown that the L^- and H^- error estimates known 
from the corresponding Euclidean setting carry over to this case. In [7] this idea is 
extended to a semi-discrete approximation of linear parabolic equations which are de¬ 
fined on a (moving) hypersurface for which the motion is a priori given by a smooth 
one parameter family of diffeomorphisms of a fixed initial surface. Here, one has to 
take care of the fact that the time derivative is defined suitably which is tackled by the 
concept of the material derivative and the spatial discretization uses a moving mesh 
(the method is called ESFEM). Furthermore, and L^(i7^)-error estimates 

are shown and in [8] the L°°(L^)-estimate is improved to the optimal order of 0{h?). 
We refer to [9] for an survey of finite element methods for surface PDFs. 

We mention further contributions to this topic in the literature. In m Runge- 
Kutta methods known from the semi-linear Euclidean setting, cf. [151 IHl HZ], are 
adapted to ESFEM to obtain a fully discrete approximation of the linear parabolic 
equation in combination with the moving surface. See also m for a backwards 
difference time discretization of this problem. In m an additional tangential motion 
of the grid for ESFEM is introduced to improve the mesh quality, or more precisely, to 
compensate a motion related possibly deterioration of the mesh. In [22] finite element 
spaces that are induced by triangulations of an ’outer’ domain are used to discretize 
partial differential equations on a surface, see also [20]. In [23] El] an Eulerian finite 
element method for solving linear parabolic partial differential equations is presented 
and a stabilized finite element method for linear parabolic equations on surfaces is 
studied. In [10] a /i-narrow band finite element method for linear elliptic equations on 
implicit surfaces is studied. See also [2] for variational problems and partial differential 
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equations on implicit surfaces. In [1] an analysis of the discontinuous Galerkin method 
for linear elliptic problems on surfaces is carried out. 

In our paper we show that the well-known L°°-estimate for the finite element 
approximation of linear elliptic equations in a two dimensional Euclidean setting, cf. 
[24] and compare also [nidads], carries over to the case of a linear elliptic equation 
on a surface in which seems to be omitted in the literature until now according to 
the author’s knowledge. 

We refer to [I] where an icosahedral discretization of the two-sphere is used to 
solve the Laplace-Beltrami equation on the two-sphere. There it is claimed (without 
detailed justification) that the quadratic order of the T““-interpolation error imme¬ 
diately carries over to an ^““-estimate of quadratic order for the discretization error, 
see the passage following Table 2 on page 1114 in [I], which is wrong, of course. 

Our paper is organized as follows. In the remaining part of the introduction we 
present the general setting and formulate our partial differential equation. In Section 
l^we introduce our notation, state some basic facts which will be used in the sequel 
several times and present the discretization of the equation. In Section |3] we present 
for completeness in the surface case the known H^- and L^-error estimates, cf. |6]. 
In Section |4| we state in Theorem 14 .1 1 our main result about the L°“-error estimate in 
the surface case and present a proof. 

We sketch the idea of the proof. We prove the estimate pointwise by using an 
approximative Green’s function v on the surface. The latter function is obtained by 
lifting a cutted-off Euclidean Green’s function from the tangent plane to the surface 
at which the appearing discrepancy to an exact Green’s function on the surface is - 
in case of relevance - supressed by the L^-error estimates from Section |31 We define 
a finite element approximation of v for which we prove an error estimate in the W^'^- 
norm which has the same order as in the corresponding Euclidean case. In doing so 
we adapt the argumentation from |24j . 

Let S' be a smooth two-dimensional, embedded, orientable, closed hypersurface 
in We triangulate the surface by a family Th of flat triangles with corners (i.e. 
nodes) lying on S. We denote the surface of class (7°’^ given by the union of the 
triangles t € Th hy Sh', the union of the corresponding nodes is denoted by Nh- Here, 
h > 0 denotes a discretization parameter which is related to the triangulation in the 
following way. Eor t G T we define the diameter p{t) of the smallest disc containing 
T, the diameter ct(t) of the largest disc contained in r and 

(1.2) h = maxp(r), 7 ^ = min [ ^ . 

reTh reTh h 

We assume that the family {Th)h>o is quasi-uniform, i.e. 7 /t > 70 > 0. We let 

(1.3) Vh = {v G C°{Sh) : v\t linear for all t G Th} 

be the space of continuous piecewise linear finite elements. 

We assume / G L^{S) and our goal is to prove error estimates for a finite element 
approximation of the unique solution u G H^{S) of the PDE 

(1.4) - Asu + u = f 

where As is the Laplace-Beltrami operator on S. In Section 0| we will assume that / 
is in addition so that u G W'^'°°{S). 

Remark 1.1. After submitting a first version of the present article to arXiv the 
author became aware of the fact that Theorem ED has been proved in m- 
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2. Notations, elementary observations and discrete formulation. Let N 

be a tubular neighborhood of S in which the Euclidean metric of N can be written 
in the coordinates (x^,x) = (x^,x’‘) of the tubular neighborhood as 

( 2 . 1 ) ga/3 = {dx°)‘^ + crij{x)dx'‘dx^. 

Here, denotes the globally (in N) defined signed distance to S and x = (a;*)i=i ,2 
local coordinates for S. 

For small h we can write Sh as graph (with respect to the coordinates of the 
tubular neighborhood) over S', i.e. 

(2.2) Sh = graph'0 = {(a;°,x) : x° = i’{x),x G S} 
where ip = iph ^ C°’^(S) suitable. Note, that 

(2.3) \D4’\<t < ch, \ip\ < ch?. 

The induced metric of Sh is given by 

(2.4) gij{ip[x),x) = 

Hence we have for the metrics, their inverses and their determinants 

(2.5) gij = aij + 0{h^), g’’^ = cr*-^ + 0{h?) and g = a + 0{h?)\aija'^^\i 

where we use summation convention. 

Let / G VE^’^(S), g G IE^’^*(S), 1 < p < oo and p* Holder conjugate of p. We 
define the so-called lift / of / to Sh by f[x) = f{ip{x), x), x G S, and correspondingly 
for g (more generally, we can do this procedure whenever we have two graphs in the 
same coordinate system and denote it by the terminus lift, furthermore, this terminus 
can be obviously extended to subsets). In local coordinates x = (x®) of S hold 

(2.6) J^{Df,Dg) = ^^cr*-^(x)v^cr(x)dx*dx^ 

(2.7) (^Df,Dg'j = ^^5*'’('0(x),x)Vg(0(x),x)dx*dx^ 

(2.8) [ {Df,Dg)= [ (Df,Dg)+0{h^)\\f\\w^,.^s)\\9\\w^..‘iS), 

Js JSh ' ' 

and similarly, 

(2.9) [ f= f f + 0{h^)\\fhiis) 

Js JSh 

where now / G L^(S) is sufficient. 

The bracket (m, v) denotes here the scalar product of two tangent vectors u, v (or 
their covariant counterparts), and later also the application of a distribution u to a 
test function v\ which meaning is on hand will be clear from context. || • ||rv'='P denotes 
the usual Sobolev norm, | • \y/k,p = X]|a|=fc ' lU'’ ■ 

Remark 2.1. Let Zq G S and TzqS be the tangent plane of S in zq. A ball 
B 2 riizo) C S, ri > 0 suitable, and a corresponding portion of Sh can be written 
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as a graph over a corresponding subset U of in a (perpendicular) Euclidean 

coordinate system ,x'^). Here, x^,x‘^ denote Euclidean coordinates in T^gS 

with center in zq and x^ denotes the coordinate axis perpendicular to T^gS so that 
TzgS = = 0}. One can consider lifts of functions between these three (pieces of) 

surfaces with respect to this Euclidean representation as graph. Analogous estimates 
to \2.!]\) . hold except for lifts from or to U , for these 0{h?) has to be replaced 

by 0(max(diam(supp/),diain(supp(7))^). Since S is compact all constants in the 
estimates and ri can be chosen independently from zq. 

Since the properties and aspects needed to prove a priori error estimates for 
finite element approximations are formulated in terms of integrals these observations 
concerning the transformation behavior of integrals essentially imply that the known 
error estimates from the Euclidean setting carry over to the surface case as far as 
convergence of at most quadratic order is concerned. 

We define 

(2.10) a : W^’P{S) x W^’P* (S') ^ R, a{u, v) = ( {Du, Dv) + uv 

Js 

and 

( 2 . 11 ) Oh-W^’PiSh) xW^’P {Sh) a{uh,Vh)=[ {Duh, Dvh) + UhVh- 

JSh 

The finite element approximation of u in (HH) is defined as the unique Uh € 14 with 

( 2 . 12 ) ah{uh,g>h)= / fh^h '^Ph&Vh 

JSh 

where fh is the lift of / to St- Note, that existence follows from uniqueness. 

Constants that do not depend on h or zq are denoted by c or co,ci if several 
appear and should be specifiable. 

Functions labelled by capital letters denote (without explicit declaration) the lift 
via the representation as graph with respect to the tubular neighborhood iV of S of 
the function denoted with the corresponding small letter to the other surface, e.g. 
assume w is defined on S then W denotes the lift of w to Sh and vice versa, i.e. if iv 
is defined on S/i then W denotes the lift of w to 5'. 

3. The iJ^-estimate and L^-estimate. For completeness we give in this Sec¬ 
tion a proof of the well-known estimates stated in Lemma 13.11 and Lemma 13.21 cf. 
0 - 

Lemma 3.1. We have 

(3-1) ||u- Uh\\m{s) < 0(li)||/||i2(5). 


Proof. Let iph € 14 arbitrary then 

||u - UhUh^s) =“(“ -Uh,u- Uh) 

=a{u -Uh,u- ^h) + a{u - Uh, ^h - Uh) 

<||u - Uh\\Hps)\\u - -Uh,^h- Uh). 
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(3.2) 


We rewrite 


(3.3) 


a{u-Uh,^h-Uh) = [ fi^h-Uh)- [ fh{^h-Uh) 

Js JSh 

.2\ 


+ 0{h^)\\Uh\\m(s)\\^h - Uh\\m(s) 
=0{h^)i\\U4HHS) + Wfh^mh - UhWh^^s) 

so that we obtain from (13.21) 

(3.4) \\u-Uh\\H^(s) < 2max(mi,m2) 
where 

(3.5) 

and 


mi = inf 

Vh£Vh 


(3.6) m 2 = 0{h )||/||l 2 ( 5 )( inf ||$;, - u\\hi(^s) + 11^ - Uh\\m{S)) 

(fh^Vh 

Let u be the lift of u to Sh then 

inf llu - $;,||i/i(s) <||u - u||iii(s) + inf ||m - iph\\m{SH) + 0{h'^)\\f\\L^s) 


(3.7) 


Vh&Vh 


Vh^Vh 


< 0 (MII/IIl^(S)- 


Putting these estimates together yields the claim. □ 
The estimate in the L^-norm can be improved. 
Lemma 3.2. We have 


(3.8) 


\\u-Uh\W(S)<0{h^)\\f\\L^(S)- 


Proof. Let w € H'^{S) be the unique solution of —Asw + w = u — Uh and Wh € 14 
the corresponding unique hnite element solution to the right-hand side u — Uh- Then 
we have 


/ {u-Uhf =a{u-Uh,w) 
JS 


(3.9) 


=a{u -Uh,w - Wh) + a{u - Uh, Wh) 

<||m - - Wh\\H^(S) 

+ 0{h'^){\\Uh\\Hps) + \\f\\L^iS))\\Wh\\Hps) 
<ch\\u - 17 ;^||/ i 1 ( 5 )||u - Uh\\L^{S) 

+ Oih^)i\\Uh\\Hps) + \\f\\L^iS))\\u - UhWmsy 


4. The L°°-estimate. We assume that / £ L‘^{S) is in addition so that u £ 
W'^’°°{S). The following theorem states our main result. 

Theorem 4.1. There holds 


I|m - Uh\\L°°{S) < log/l|||M||w2.“(S)- 


(4.1) 
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The proof of the corresponding Euclidean statement is well-known, cf. 

The purpose of the remaining part of this section is to prove Theorem 14.II 
Let zq G S and 

(4.2) (ps-.U ^ B 2 r^izo), {0,x) ^ x = {x\x‘^) 

be the representation as graph of B^r^izo) C S over U C according to Remark 

[211 

Definition 4.2. We set Bj = Bjrj^(zo), j = 1,2, p = ps v be the 

Euclidean Green’s function with respect to —A + I in T;,,,S = with singularity in 
zq, i.e., more precisely, 

-Av + v=dzo in Bioo{zo) C 
dnV=0 ondBioo{zo). 

Let ( G C^{B3 ^_^(zo)), f = 1 in Bi, be a cut-off function and set 

(4.4) u(a;) = v((/3"^(a;))C(a:), l{x) = l{p~^{x))Cix), x G B 2 , 

where l{z) = log \z — Zo\. 

There holds 

(4.5) y-^i^H^^Bioo{zo)), 

ZTT 

cf. [24l Lemma 1, page 687]. W.l.o.g. we may assume ||t||c2 (B2 \Bi) — and 

(4.6) Wi-^l\\w2.2(^s) <c 
where co,c are independent from Zq. 

The next Lemma shows in which sense v is an approximative Green’s function. 
Lemma 4.3. Let v be as in Definition \4.2\ Let w G L[^{S) fl C^{S) then 

(4.7) I {-Asv -\-v,w) - w{zo)\ < c||u;||l 2 ( 5 ). 


Proof Let rj G C^{Br,{zo)), 1713 , 1 / 2 (^ 0 ) = 1^1 ^ ^ and \Dp\ < c. We write 


(4.8) 


= r]w -I- (1 — ri)w = ' 


and have (in local coordinates induced from p) 


(4.9) 


{—Asv -j- v,w) = (^—Asv V, w^') {^—Asv -I- v, w'^') 




/ ■ dv dw"^ ~ 2 \ r- 

» + ™ 1^- 


We rewrite the second integral on the right-hand side of (14.91) as 


92 / 


~v . ^ 2 f d i, dv 2 


(4.10) 


dv 1 


d 


9- 


dx’- 2 

=0(1)||u;||l2(s) 




d'^v 

dx^dxl 
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where we used integration by parts and all integrals are over U \(p (^o))- 

Let r denote the distance to Zq in T^qS then we obtain 


\S'^ -g'i\ <cr^, 

\g^^\<c 


d 

dx^ 

(4.11) 

dl 

C 

d^l 


dx' 

< 

r 

dx'dxi 


We rewrite the first integral on the right-hand side of (14.91) as 


(4.12) 


dv dw^ 
dx'- dx^ 
dv dw^ 

a"-' —r-r 

dx'^ dx^ 
ijj \ dx'- dxi 


+ VW ] ^ 


+ VW {^/g-l) + / (ff* 


• VW 


=w(2;o) + 0(l)||w||i2(s) 




dv dw^ 
dx' dx^ 


where we used (j4.1ip . Holder’s inequality, ()4.3I1 and that we are allowed to perform 
integration by parts in the integrals with the factors {y/g — 1) and g'^ — 5'^, note, that 


(4.13) 





d^-u 

dx'dx^ ^ 


< 


_ g^J) 


<c|k||L 2 (S)- 


d'^{v-l) 

dx'dx^ 


d^l A 1 

dx'dx^ y ^ 


□ 

Remark 4.4. From now, we denote the approximative Green’s function v by g 
(there will be no ambiguity with the symbol for the determinant of the metric). 

We define an approximation g^ G 14 of g by 

(4.14) ah{gh, Vh) = a{g, 14) 

for all Vh GVh- 

Lemma 4.5. Assume 

(4.15) \\g - GhWw'-AS) < ch\\ogh\ 

then Theorem \).l\ follows. 

Proof. From Lemma 14.31 we conclude 

{u-Uh){zo) =a{g,u - Uh) + 0{h^)\\f\\L2f^s) 

(4-16) =ah{G — gh,U — Uh) + 0{h'^)\\g\\w-^v(s)\\u — Uh\\w'’’=°{S) 

+ ah{gh,U - Uh). 

We estimate 

ah{gh, U — Uh) =a{Gh,u) — ah{gh,Uh) + 0(/i^)|lG?i||vi/i,i(s)||u||v^i.oo(s) 

=0(/l^) (||Gh||wi.i(S)||u||vyi.oo(s) -I- ||G?i|4i(s)||u||v^2.oo(s)) . 
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(4.17) 























Furthermore, we have 


(4.18) ah{G - gh,U - Uh) =ah{G - gh,U - Vh) + ah{G - gh,Vh - Uh), 
rewrite the second summand by using Lemma 14.31 as 
^h{G gfn Vh Uh') 

= 0 '{g, Vh - Uh) + 0{h'^)\\g\\w^,i{s)\\Vh - Uh\\w^,^(s) 

^h\9hi '^h '^h) 

<0(h^)||g||wi'i(s)||14 - Uh\\w^’°°{S) 
and estimate the first summand as follows 

(4.20) \ah{G - gh,U - Vh)\ <||ff - G/i||wi,i(s)||u - 14||vi/i,oo(5) 

We let Vh S Vh be the interpolation of u and obtain the claim from 

I|m - [//tllwi.oo(s) <I|m - Vh\\w^.<^is) + \\Vh - Uh\\w^’°°{S) 

(4.21) <||m - 

+ ch ^ (\\Vh — u||wi. 2 (s) + \\u — UhWw^’^iS)) 

which holds in view of estimate ()4.37p . □ 

In order to show ()4.15p we prove as first step the following Lemma. 
Lemma 4.6. We have 


(4.22) \\9-Gh\\LHS)<ch 
where c is independent of Zq. 

Proof. Let r be a triangle in Th containing the lift zq of zq, and let g be the linear 
function with 

(4.23) J qp=p{zo) 

for all linear functions p. Because r contains a disk of radius 70/1 we see that 

(4.24) sup|( 7 | < c/i-2. 

T 

We extend the domain of definition of q to Sh by zero and set S = Q. 

We define 


(4.25) ^P{v)=a{g,v)-{S,„v), v G W^’^{S) n G°{S), 

From Lemma 14.31 we deduce that 


(4.26) |■^/'(^^)| < co||u||l 2 ( 5 ) 

so that by Hahn-Banach Theorem ip can be extended to a linear functional on Lf{S)- 
denoted by ip as well - with norm < cq. W.l.o.g we may assume that ip G Lf{S). 

Let w G H^{S) and Wh G Vh be the solutions of 


(4.27) 



Vv G ff^(S) 


ah{wh,Vh)= 'i'Vh 'ivhGVh- 
JSh 
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Lemma o leads to 
(4.28) 


\\w - WhWmiS) < ch 


where c independent from Zq. 
Let Zfi € Vh with 


(4.29) ah{zh,Vh) = {Sio,Vh) =Vh{zo)- 
Let g solve 

(4.30) -Asg + g = 5. 

Since Zh can be seen as finite element approximation of g we have in view of Lemma 

o 

(4.31) \\g-Zh\\H‘(S) <ch'^~"\\HL'^(S) s = 0,1. 

In view of 


ah{gh-Wh,Vh) =a{g,Vh) - j 'i’Vh 

JSh 

=a{g.Vh)- [ Wh + 0{h^)mL.^s)\\V4LHS) 

Js 

= {S.„Vh) + 0{h^)\mL2^s)\\Vh\\LHS), yvnGVh 

we deduce 


(4.33) \\zh - {gh - wh)\\L^(s) < 0{h^). 

To estimate 


(4.34) ||g - G/i|1l2(5) = \\{g-w-g) + {g- Zh) + {w- Wh)\\L'^{s) + OQi^) 
we need to estimate \\g — w — g\\L 2 (^gy Let y) G C°°{S) and w a solution of 

(4.35) — Agw + w = (fi 


then 


(4.36) 


/ {g-w- =a{g -w-g,w) 

Js 


= (Szn - S.w) 


= {5z„-5,w-wi 


<||n; - ■W7 ||l.»(s) + Il^l|i2(s)l|w - w/||l 2 (s) 

<0{h)\\ip\\L2(^s)- 


Here, Wj denotes the linear interpolation of w, wj its lift to S and we used, cf. [31 
Theorem 4.4.20], 


(4.37) hi\\x-Xi\\L-^{S)+Yl h^Wx- Xi\\w^’P(S) < ch‘^\\x\\w^’P{S), 1 < P < oo, 

3=0 
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for X G H'^{S) and xi the linear interpolation of x (, and the right-hand side possibly 
unbounded). □ 

Remark 4.7. (i) Estimate j-j-lS]} follows immediately if we show 

(4.38) ll/'-L,,|| ivi,i( 5 ) < c/i|logh| 
where Ih G 14 is defined by 

(4.39) ah{lh,Vh) = a(i, Vu) 'ivu G Vu- 
(ii) There holds 

(4.40) Ilf- L/j||i 2 ( 5 ) < ch. 


Proof, (i) We have 

(4.41) Qhigh - -^h,Vh) = a{g - ^4) 

so that in view of (j4.5|) we conclude from Lemma [34l that 

(4-42) \\Gh - - {g- ^OIIri(S) < ch 

and the triangle inequality implies (j4.15l) . 

(ii) Use (I4.42|) . the triangle inequality and Lemma [4.61 □ 

In the remaining part of this section we prove (I4.38|) . We recall that l{z) = 
log \ z — zq\ is defined in Tz„S = that r denotes the distance to zq in T^gS and 
state that I has bounded mean oscillation in the following sense. 

Lemma 4.8. Let zi G and 0 < p < oo. Then there is a constant 4 G R 
depending on Zi and p such that 

(4.43) [ {l-lo?<97rp\ 

j {\z-Zl\<p} 

Proof. This is the assertion of [24l Lemma 2 on page 688]. □ 

Remark 4.9. In the following we will consider lifts of objects defined on i?2ri (-Zo) C 
S, U C. TzgS or a suitable portion of Sh to another of these three surfaces with respect 
to the representation as graph over U in (perpendicular) Euclidean coordinates as de¬ 
scribed in Remark \2.1\ By adding the superscripts S, T or h we indicate to which 
surface the object is lifted, e.g. let M C B 2 ri{zo) C S then M'^ denotes its lift to 
TzgS. Similar correction terms as in 112.^) and i2.9\) appear when we lift integrands 
of (with a power ofr) weighted -norms, i.e. if we estimate such a norm then the 
lift produces (at most) a constant as factor on the right-hand side of the estimates. 
We estimate the error E = I — Lh near zq. 

Lemma 4.10. Let 0 < p < cih be given and B = {\z — zo\ < p} a ball in T^gS. 
Then 

(4.44) / \pf^{z) - zof\DE\P < cp^h^-P 

JsnBS 

for 1 < p < fj -\-2. 
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Proof. We have \Dl\ < ^ where r = — zq\. In view of Remark l4.9l we may 

w.l.o.g. consider r as a function as well on and B = B^ and get 


(4.45) 


Is 


SnBS 


\^.,{z)-Zof\DE\P <c 


L 


isnBS 


r^{\Dl\P + \DLhn 


s'* 


'\Dk\- 


By Lemma 14.81 there is Zq S R so that 

(4.46) \\l — lo\\L^{B'r) < ch. 

We get (using an inverse estimate to bound a by a L^-norm) 


[ r^\Dh\P=[ r^\D{k-lo)\^ 
Jb>' Jb>' 


(4.47) 


<cp^h^ sup \D{lh - Zo)r 
s'* 

<cp^hf‘h~'^P\\lh - Zo|li2(Bh) 

<cp^h~'^P{\\Lh - ^11^2(53) + ||Z - lid\\L^(BT)Y 
^cpI^h^-P 

in view of ()4.40|) and (I4.46p . □ 

Remark 4.11. If we choose fi = 0 in Lemma\4Bd\we obtain that \\DE\\Bi(snBS) = 
Oih). 

We estimate the error E outside B^, B as in Lemma [4. 101 which means de facto 
in i ?2 \ since suppZ C B3^^{zo) and get 


(4.48) 


and 


(4.49) 


/ \DE\<([ r-^y ( / 

Jb2\BS \Jb2\BS / \Js J 


< c|log/i|5 / r‘^\DE\ 


[ r^\DE\^ = [ {DE,D{r^E)) -2Er{DE,Dr) 

Js Js 


< I {DE,Dir^E)) + 2[ j E^ 


'•\DEf 


which leads by Peter-Paul inequality to 

(4.50) / r'^\DE\^ <2 / {DE,D(r'^E)) + A / E^ <2 {DE,D{r^E)) + ch^ 

Js Js Js Js 

in view of Remark 14.71 (ii). The next goal is to show 

f {DE,Dir^E)) <1- [ r^\DE\^+ ch^\\ogh\ 

Js 4 Jg 
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(4.51) 










which implies (14.381) . Define 


(4.52) = {t & Th ■■ dist {zq,t) >h}, Di = [J r C Sh- 

tGTI 


and note, that for small h 

(4.53) {z € S : dists( 2 :, zq) > 3h} C Df 

and I G Let Ij be a function in 14 which equals I at all nodes in Di. Let Ij 

denote the lift of // to S. 

Lemma 4.12. There hold 

(4.54) [ [Ti-lf<ch\ [ r-^\D{r^{Ii-!))f <ch^\\ogh\. 

Jnf jQf 


(4.55) 


Proof. Let t G then 

\\l ^ch ||Z||^y2,oo(.rS) 

<c/i^“®(minr)“^ s = 0,1. 


Since min,- r > h, maxi- r — min,- r < h we have 

max,- r 


(4.56) 

and hence for fi >0 


< 2 


mm,- r 


(4.57) f r^{I — liY + hf [ r^\D{l — lh)f < c [ r^h‘^{mmr) ^ < c / 
JtS JtS JtS JtS 

Summing over all t G implies the Lemma since 


(4.58) 

and 

(4.59) 

□ 


[ r^-^h^< 
JtS 


Jch^+^, if/3<2, 

1 c| \ogh\h‘^, if /3 = 2 


r-^\D{r^{l - 4 ))r < 8 (Z - hY + 2r"|D(Z - li)\f 


We conclude 


(4.60) 


h - lh\\L°^{Q.i) <ch \\h - lh\\L2{Qi) 

<ch ^{\\Ti — ^||l 2 (qS) + \\Lh — ? 1 |l 2 (qS)) 
<c 


in view of Lemma [4.121 and Remark l4.7l fiil. 

Lemma 4.13. Let € 14 and v = (r^tp)/ G 14 the linear interpolation ofr^ip in 
Di then 

(4.61) / r-^\D[r^ip-v)\^ <c ( 
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Proof. For t £ we have 


(4.62) 

because D‘^{(f\T) 

(4.63) 


\r'^(p-v\w^,^(r) <ch\r'^(p\w2.<=-(r) 

2 

<ch E 

i=i 

= 0. In view of (14.561) and r > /i on fli there holds 
|r^|w3,°=(r) < cinfr^”-^ < cinirh^~^ 


and in view of an inverse estimate 


(4.64) \^p\w2-j,^(^r) <ch^ ^WpWl^t)- 
Applying these estimates in (14.621) gives 

(4.65) \r‘^<P - <c/i"^ inf r||(/3||i2(^) 

which leads to 


(4.66) 


J r ^|ZI(r^i^ — u)|^ < c J 


by estimating the integrand in the L°°-norm. Summing over t G gives the claim. 


Lemma 4.14. Estimate holds. 

Proof. For Vh G 14 we have 


(4.67) 


a(E,Vh) = 0{h^)\\Lh\\H^iS)\\yh\\H^{S) 


and estimate 


f {DE,D{r^E))= f {DE,Dir^E-Vh))- j EVh 
J s ^ s ^ s 

+ 0{h^)\\Lh\\m{s)\\yh\\m{s) 

r 

/ {DE, D{r^E - Vh)) + c{h^ + h\Vh\w^.^^s\n,)) 

«/ r2i 


(4.68) Lemma \ 4 :. 10 \ Remark 14.71 
< 


+ / yh + Oih‘^)\\Lh\\m{s)\\yh\\m{s)- 

Js 
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If Vh interpolates r‘^{li — Lh) in f2i then 


(4.69) 


f {DE,D{r^E-Vh)) <^ [ r^\DE\^ +4 f r-^\D{r^ E - Vh)\^ 
inf 16 is inf 

[ r^\DE\^+8 f r-^\D{r^{l-Ii))\^ 

16 is inf 

+ 8 [ r-^\Dir^Ii-Lh)-Vh)\^ 

inf 

Lemma 14.12^2). Lemma 14.131 1 /* /* 

^ / r2|i4i5|2 + c/i2|log/i|+c / (rz-L^)2 
16 is inf 


< 


-i^ / I+ c/i^I log/iI 
16 is 


+ c / ih — + {I — LfiY 

\inf inf 


Lemma 14. 12| i^emarfc EZli 
< — 

16 ./s 


f r^\DE\^ + ch^\\ogh\. 
is 


We use (I4.69|) to estimate the first summand on the right-hand side of (j4.68|) and 
obtain 


(4 70) Is Tq Jg + c/i| 14 |^yi,oo( 5 \QS) -I- 

+ 0{h'^)\\Lh\\m(s)\\Vh\\mis) 

We estimate Vh with standard interpolation estimates 


h\vh\w^,^{Sh\ni)+h <c sup \r'^{li-ih)\ 

s\nf 

(4.71) =sup|r^(r/- L?,)| 

snf 

<ch? 


where we assume w.l.o.g. that Vh is zero at all nodes in the interior of Sh \ fli and 
for the last inequality estimate (14.601) . Furthermore, we have 

lk?i||L2(ni) <c\\r‘^(ii - Lh)\\L2^QS) 

<c\\Il — ?||s2(f2f) + \\l — -b/l|lL2(QS) 

<ch 

in view of Lemma 14.121 and Remark 14.71 and 

(4.73) I114||//i(s) < ch~^\\Vh\\L^(s) < ch 
and 

(4.74) l4h||i/i(S) ^ ch~^\\Lh\\L^{s) < c. 

□ 
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